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SUMMARY 


Numerical solutions are presented for three-dimensional laminar and turbulent 
flow in curved ducts of rectangular cross section and significant curvature. The 
analysis is based on a primary-secondary velocity decomposition in a given coordi- 
nate system, and leads to approximate governing equations which correct an ^ priori 
inviscid solution for viscous effects, secondary flows, total pressure distortion, 
heat transfer, and internal flow blockage and losses. Solution of the correction 
equations is accomplished as an initial-value problem in space using an implicit 
forward-marching technique. The overall solution procedure requires significantly 
less computational effort than Navier-Stokes algorithms. The present solution 
procedure is effective even with the extreme local mesh resolution which is neces- 
sary to resolve near-wall sublayer regions in turbulent flow calculations. Computed 
solutions for both laminar and turbulent flow compare very favorably with available 
analytical and experimental results. The overall method appears very promising as 
an economical procedure for making detailed predictions of viscous primary and 
secondary flows in highly curved passages-. 
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INTRODUCTION 


An economical method for predicting three-dimensional turbulent subsonic flow 
in curved flow passages is of considerable value in turbomachinery applications. Of 
particular interest here is the flow and heat transfer in curved passages including 
those shaped like the flow passage between adjacent blades of a turbine. In this 
application, the ultimate goal is the detailed prediction of secondary flows, heat 
transfer and aerodynamic losses in turbine blade passages. This problem is of par- 
ticular importance for gas turbine design, since endwall cooling becomes a critical 
problem as turbine inlet temperatures are increased to achieve high cycle efficien- 
cies. It is generally believed that heat transfer to the turbine endwall and endwall 
losses are strongly influenced by the large secondary flows which exist in turbine 
airfoil passages. These secondary flows, which are caused by the turning of the 
primary flow, can remove insulating endwall boundary layer fluid and replace it with 
hot mainstream gases, thereby increasing endwall heat transfer. The secondary flows 
also represent angular momentum not recoverable as thrust and thus are losses which 
are of prime concern to the turbine designer. 

A definitive approach for computation of three-dimensional flow in curved pas- 
sages would be numerical solution of the Navier-Stokes equations. Although feasible, 
in three dimensions this approach requires considerable computational effort and is 
most attractive when no viable alternative exists. The present study considers an 
alternative approach for predicting such flows and is a development of the primary- 
secondary velocity decomposition method of Briley and McDonald [1] for application to 
viscous subsonic flow in smoothly curved geometries. The objective of this approach 
is to introduce approximations which adequately represent essential physical processes 
of interest and yet lead to governing equations which can be solved much more econom- 
ically than the Navier-Stokes equations. In the present treatment of developing 
flows in curved passages, it is presumed necessary to provide an adequate representa- 
tion of primary flows, secondary flows, viscous effects, and their local interactions. 
In the approach taken here, an inviscid flow approximation is first obtained for the 
geometry in question. The inviscid flow may, for example, satisfy an elliptic govern- 
ing equation requiring downstream boundary conditions (e.g. a velocity potential 
equation) and thus include transverse variations in streamwise pressure gradient 
usually associated with flow in curved passages. This a. priori inviscid flow is 
then corrected by solution of a set of approximate governing equations and boundary 
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conditions which constitute a well-posed initial-value problem in space. These 
correction equations apply in both viscous and inviscid regions and account for 
secondary flows, viscous effects, heat transfer, total pressure distortion, and 
internal-flow blockage and losses. Solution of the correction equations as an 
initial-value problem provides the desired reduction in computational effort of one 
or more orders of magnitude over Navier-Stokes solution procedures. 

The present development of this approach includes an improved solution algorithm 
which provides for solution of a pair of correction equations (governing vorticity 
and vector potential) as a coupled system, using an iterative linearized block 
implicit (LBI) scheme. This improvement permits implicit specification of no-slip 
boundary conditions for the secondary flow velocities and removes a previous assump- 
tion that flow in the near wall region is collateral. The method is used to compute 
laminar and turbulent flow in curved ducts of rectangular cross section, and these 
results are compared with experimental measurements and other numerical predictions. 
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RELATED WORK 


Secondary Flow Theory 

The concept of dividing the flow into primary and secondary components is 
central to the present formulation of an initial-value correction to an ^ priori 
inviscid flow approximation. The presence of strong secondary flows is also of 
particular Importance in applications of interest here, and the insight afforded by 
secondary flow theory is thus of considerable value. Excellent reviews of secondary 
flow theory and its applications are given by Hawthorne [2,3] and Horlock and 
Lakshminarayana [A]. A key and relevant feature of many studies based on secondary 
flow theory is that they provided a basis for treatment of strong rotational but 
inviscid secondary flows as an initial-value problem, without the necessity of a 
three-dimensional iteration procedure (cf. Rowe [5]). Previous approaches based on 
secondary flow theory have generally neglected viscous effects, however, and instead 
have introduced a somewhat arbitrary "cut off" velocity to exclude the viscous 
region. Hawthorne [2] has pointed out that ambiguity of the "cut off" velocity can 
have a significant effect on the predictions obtainable from the inviscid theory. 

The present approach includes a detailed representation of viscous effects on both 
primary and secondary flows, and since no distinction between inviscid and viscous 
flow regions is required, there is no need to introduce a "cut off" velocity. This 
is accomplished by introducing approximations into a generalized (exact) expression 
from viscous secondary flow theory given by Lakshminarayana and Horlock [6]. In 
addition to viscous effects, the present treatment also accounts for generation of 
secondary vorticity by inviscid mechanisms such as turning of streamlines, for trans- 
port of secondary vorticity, and for distortion of the primary flow. 

Forward-Marching Algorithms 

In recent years, a variety of numerical methods based on forward marching al- 
gorithms have been devised for approximating three-dimensional viscous flows for 
which one "primary" velocity component does not change sign (separate) . These 
methods are often characterized by terminology such as "parabolized Navier-Stokes" , 
"parabolic flow", "partially-parabolic flow", etc., because elements of either the 
solution procedure or some of the equations solved or both have "parabolic" or 
"elliptic" qualities. Unfortunately, this terminology does not identify the relevant 
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physical approximations made or their significance, and does not distinguish these 
approximations from properties of the solution algorithm and the differential or 
difference equations. Different methods within the same "category" will in some 
instances give significantly different results. 

Here, methods are identified as either initial-value methods or methods which 
require three-dimensional iteration, and an attempt is made to identify the most 
important assumptions made. Although some type of forward-marching algorithm may be 
used in both the initial value and global iteration approaches, the present distinc- 
tion is intended to emphasize that initial-value methods provide a solution after a 
single forward-marching integration and normally require significantly less computa- 
tional effort than global iteration methods. Global iteration methods on the other 
hand provide a solution to the system of governing equations only upon convergence 
of the iterative process. This distinction between the two categories is important, 
since the approximations made in a proposed method should permit a reduction in com- 
putational effort which justifies their use vis -a- vis solving the Navier-Stokes 
equations without approximation. As a comparison familiar to the present authors, 
recently developed split linearized block implicit (LBI) algorithms [7] have achieved 
CDC 7600 run times on the order of 5 minutes per 1000 grid points for the compressible 
Navier-Stokes equations in three dimensions [8]. The initial-value correction method 
employed in the present study requires about 20 seconds per 1000 grid points, a 15:1 
savings . 

Initial Value Methods 


Patankar and Spalding [9] and Caretto, Curr, and Spalding [10] formulated 
equivalent sets of approximate governing equations for use in treating a class of 
three-dimensional viscous flows without primary flow reversal. They derived these 
governing equations by neglecting streamwise diffusion and assuming that streamwise 
pressure gradients are independent of the transverse coordinates. This assumption 
can be viewed as equivalent to adding a mean correction to pressure gradients from 
an inviscid flow with constant velocity, and is thus appropriate for flow geometries 
without significant curvature. In Refs. [9] and [10], numerical methods for treating 
these equations were proposed. The numerical method of Patankar and Spalding [9] has 
since been widely used and introduces a further assumption that streamwise increments 
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in computed transverse pressure gradients are proportional to transverse-velocity 
corrections introduced to satisfy continuity. The particular pressure-velocity 
relationship is obtained by neglecting all off-diagonal terms in finite difference 
formulas which approximate the transverse momentum equations. As a consequence, the 
computed transverse pressure distribution does not satisfy a (numerically) consistent 
difference approximation of the transverse momentum equations evaluated with con- 
tinuity-corrected velocities, except for fully-developed flows for which velocity 
and pressure corrections are uniformly zero. 

As a means of treating flow geometries with significant curvature, Briley [11] 
suggested an alternative formulation in which streamwise pressure gradients are ap- 
proximated as the sum of gradients from a potential flow solution in the curved 
geometry and a mean correction which is independent of transverse coordinates. Here, 
the assumption is that the potential flow is to be corrected and that the mean 
pressure drop computed as part of the initial value correction process does not vary 
across the duct cross-section. In the numerical method of Briley [11] , it is assumed 
that transverse-velocity corrections introduced to satisfy continuity are irrotational ; 
however, the transverse momentum equations are re-solved to obtain the transverse 
pressure distribution consistent with the continuity-corrected velocities. 

Ghia and Sokhey [12] have extended the formulation of Patankar and Spalding [9] 
and numerical method of Briley [11] to treat incompressible flow in straight or cir- 
cular arc ducts with either rectangular or polar cross sections. They were successful 
in obtaining solutions for laminar flow development in circular arc ducts of rectang- 
ular cross section and strong curvature. 

As part of a comprehensive study of developing flow in straight rectangular 
ducts, Rubin and Khosla [13] derived from asymptotic theory a system of governing 
equations based on velocity, vorticity, and a potential and stream function for the 
crossflow. These equations are valid in the "fully viscous" region also addressed by 
the initial-value methods of [9-12]. The Rubin-Khosla governing equations are in- 
dependent of Reynolds number, and are solved numerically as an initial value problem 
without further approximation. 

The present study is a development of the initial-value correction approach of 
Briley and McDonald [1] and provides a general method for treating flows in smoothly 
curved geometries which may have significant streamwise curvature and turning. This 
method is based on a primary-secondary velocity decomposition in a given coordinate 
system. Unlike most previous methods, the primary flow velocity component may or may 
not be aligned with the marching coordinate direction. The primary flow is governed 
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by a momentum equation in which streamwise pressure gradients are approximated from 
an a priori inviscid-flow analysis, computation, or approximation. The secondary 
flow is uniquely determined, through scalar and vector surface potential calculations, 
from the primary velocity and secondary or streamwise vorticity, the latter of which 
is computed here from an approximation to an exact equation governing viscous secon- 
dary flows. This system of equations corrects an a priori inviscid solution for 
viscous effects, secondary flows, total pressure distortion, heat transfer and in- 
ternal flow blockage. Solution of the correction equations is accomplished as an 
initial-value spatial marching process. The dependent variables (primary velocity, 
secondary vorticity, scalar and vector surface potentials) are analogous to those 
used by Rubin and Khosla [13] in Cartesian coordinates. 

Other developments of this particular initial value correction method have been 
applied to other problems. Shamroth and Briley [14] have obtained predictions for 
the airfoil tip vortex generation process. Kreskovsky, Briley and McDonald [15] have 
obtained predictions for turbofan exhaust lobe mixer flows, which requires treatment 
of turbulent mixing in the curved interface of two flow streams having large dif- 
ferences in stagnation temperature and flow angle, and slight differences in stagna- 
tion pressure. Further results for this problem have been computed by Anderson, 
Povinelli, and Gerstenmaier [16]. The analysis has been reformulated for use with 
constructed coordinate systems by Levy, McDonald, Briley, and Kreskovsky [17] and 
applied to duct and diffuser geometries having curved centerlines and superelliptic 
cross sections. 

Global Iteration Methods 


Methods within this category neglect streamwise diffusion and employ forward- 
marching algorithms as part of a global iteration process. Although streamwise dif- 
fusion is neglected, the system of governing equations being solved remains elliptic 
(for subsonic flow) and is thus subject to one or more downstream boundary conditions. 
To permit the use of forward-marching algorithms under these circumstances, an im- 
balance is introduced in the system of governing equations during the forward-marching 
process, either in continuity or in the representation of pressure gradients. Cor- 
rections are introduced after each iteration to reduce this imbalance, and the system 
of equations and boundary conditions is satisfied once the global iteration converges. 
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The "partially parabolic" method of Pratap and Spalding [18] is a global itera- 
tion which begins with a "guessed" pressure field and performs iterated forward 
marching sweeps of the flow field. During these forward-marching iterations, the 
three momentum equations are solved for velocity, and two-dimensional velocity and 
pressure corrections are determined to remove a local imbalance in continuity. 

Various strategies are then employed to modify the three-dimensional pressure field, 
so as to improve the local continuity balance on subsequent iterations. Moore and 
Moore [19] have proposed a somewhat related global iteration method which, given an 
estimated pressure distribution, computes velocities and a local two-dimensional 
pressure correction by a forward-marching algorithm. After each forward marching 
iteration, a three-dimensional global pressure correction is computed from this local 
two-dimensional pressure correction for use in the next iteration. To provide im- 
proved understanding of the predictions and possibly increase the rate of convergence, 
Moore and Moore [19] suggest converging the global iteration in three separate stages, 
each of which is itself a global iteration. The complexity of the flow model is 
increased in each stage; first is inviscid flow with uniform inflow, second is 
inviscid flow with nonuniform stagnation pressure at inflow, and third is viscous 
flow. Upon convergence such that all pressure and velocity corrections are uniformly 
zero, it appears that the only assumptions made by the global iteration methods of 
Pratap and Spalding [18] and Moore and Moore [19] are the neglect of strearawise dif- 
fusion, the condition of no reversed flow, and approximations introduced through 
boundary conditions and turbulence modeling. 

Dodge [20] has suggested a global iteration method which divides the velocity 
vector into rotational and irrotational portions which become defined only upon con- 
vergence of the global iteration and upon selection of a set of boundary conditions 
governing the rotational-irrotational velocity division. A key assumption in Dodge’s 
method is that all pressure gradients depend solely on the irrotational component of 
the velocity vector and are independent of the remaining rotational component. In 
view of this assumption, the method of Dodge does not seem suitable for applications 
having regions of rotational but inviscid secondary flow. Dwoyer [21] has also 
pointed out that the boundary conditions defining the velocity division are arbitrary 
and has shown that selection of these boundary conditions is crucial. 
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ANALYSIS 


The governing equations are derived through approximations made relative to a 
curvilinear but not necessarily orthogonal coordinate system fitted to and aligned 
with the flow geometry under consideration (cf. Fig. 1). The coordinate system is 
chosen such that the streamwise or marching coordinate either coincides with or is 
at least approximately aligned with a known inviscid primary flow direction, as deter- 
mined for example by a potential flow for the given geometry. Transverse coordinate 
surfaces must be either perpendicular or nearly perpendicular to solid walls or 
bounding surfaces, since diffusion is permitted only in these transverse coordinate 
surfaces . 

Equations governing a scalar viscous correction u_^ to a known inviscid primary 
flow velocity U ^ , and a secondary vorticity normal to transverse coordinate sur- 
faces are derived utilizing approximations which permit solution of the correction 
equations as an initial-value problem, provided reversal of the composite streamwise 
velocity does not occur. Terms representing diffusion normal to transverse coordi- 
nate surfaces are neglected. Approximate pressure gradients are derived from the 
inviscid primary flow and Imposed in the streamwise momentum equation. These pres- 
sure gradients are the sole means of accounting for the elliptic influence of down- 
stream boundary conditions in curved flow geometries. The secondary vorticity is 
computed as an approximate application of viscous secondary flow theory. Secondary 
flow velocities are deteirmined from scalar and vector surface potential calculations 
in transverse coordinate surfaces, once the primary velocity and secondary vorticity 
are known. 

Primary-Secondary Velocity Decomposition 

In what follows, vectors are denoted by an overbar, and unit vectors by a caret. 

The analysis is based on decomposition of the overall velocity vector field U into a 

primary flow velocity U and a secondary flow velocity U . The overall or composite 

p s 

velocity is determined from the superposition 

U = Up + Uj. (1) 
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The primary flow velocity is represented as 


Up 


( 2 ) 


where is a known inviscid primary flow velocity satisfying slip conditions and 
determined for example from an ^ priori potential flow solution in the geometry 
under consideration. The (non-dimensional) scalar quantity u^ is a viscous velocity 
profile factor which introduces viscous shear layers and may also correct for inter- 
nal flow blockage effects. The viscous velocity correction u^ is determined from 
solution of a primary flow momentum equation. The secondary flow velocity is 
derived from scalar and vector surface potentials denoted cp and respectively. If 
i^ denotes the unit vector normal to transverse coordinate surfaces (also presumed 
here to be in the direction of the marching coordinate), if p is density, and if 
is an arbitrary constant reference density, then is defined by 

D, 5 V^4> + (Pg/p) 


where is the surface gradient operator defined by 

Vs ^ 

It follows that i- • U lies entirely within transverse coordinate surfaces. 

Equation (3) is a general form permitting both rotational and irrotational secondary 
flows and will lead to governing equations which may be solved as an initial-boundary 
value problem. The overall velocity decomposition (1) can be written 


U = UjUy + VjC/) + VXfjl// 


(5) 


Surface Potential Equations 


Equations relating (j) and ij; with u^, p, 
can be derived using Eq . (5) as follows: 


and the secondary vorticity component 
From continuity. 


v-pu = 0 = V-pU^U^ + 


( 6 ) 
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and from the definition of the secondary vorticity, 


i|. VX U H 


= I, 


VxUjUy 


+ i ■ Vx(o /p)Vx I d, + { • VXV <#) 


(7) 


Since the last term in each of Eqs. (6,7) is zero by vector identity, Eqs. (6,7) can 
be written as 




( 8 ) 


i,-Vx(^^/p)Vxi,v/, = n, - i,-Vxu,u^ 

The last term in Eq. (9) is identically zero in a coordinate system for which and 
have the same direction, and is small if i^ and are approximately aligned. 
Further simplification is possible if is an incompressible potential flow satisfy- 
ing V • Uj = V X = 0, since the last term in each of Eqs- (8) and (9) can then be 
written - • Vpu^ and - i^ • x respectively. In any event, given a know- 

ledge of u^, f2, and p, the surface potentials (j> and can be determined by a two- 
dimensional elliptic calculation in transverse coordinate surfaces at each streamwise 
location. In turn, can be computed from Eq. (3), and the composite velocity U 
will satisfy continuity. Equations for u^ and are obtained from the equations 
governing momentum and vorticity, respectively. 


Primary Momentum and Pressure Approximation 

The streamwise momentum equation is given by 

‘,-[(U-V)0 +(Vp)//>] -?,.F 

where p is pressure and pF is force due to viscous stress. Terms in F representing 
streamwise diffusion are neglected; however, since the viscous terns are complex for 
compressible flow, the modified viscous force is temporarily denoted F', and further 
consideration of viscous terms is deferred to a later section. 
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The remaining assumption for Eq. (10) concerns the pressure gradient term and 
is designed to permit numerical solution as an initial value problem. The assump- 
tion that the streamwise pressure gradient is independent of transverse coordinates 
is inadequate for curved flow geometries, and instead the streamwise pressure gradi- 
ents are approximated from an inviscid flow obtained by an a priori analysis, com- 
putation, or approximation, and are corrected for viscous effects by a mean pressure 
gradient term which depends only on the x^ coordinate. An obvious choice for an 
inviscid velocity is the potential flow for the geometry under consideration, 
although other choices for may be used. The inviscid flow may satisfy an elliptic 
governing equation requiring downstream boundary conditions (e.g. a velocity poten- 
tial equation) and thus include transverse variations in streamwise pressure gradient 
usually associated with flow in curved passages. 

It is possible to consider the use of pressure gradients from a rotational in- 
viscid flow determined by solution of the Euler equations. It should be pointed out, 
however, that unlike the scalar potential equation, the computational labor of solv- 
ing the Euler equations is commensurate with that required for solution of the full 
Navier-Stokes equations with a similar number of mesh points. Solution of the Euler 
equations may require fewer grid points since thin shear layers are omitted and 
hence do not require resolution by the grid used. On the other hand, omitting the 
shear layers entails further approximation in selecting a "cut-off” criterion 
between inviscid and viscous regions at inflow (cf. Hawthorne [2]), and the predicted 
flow may be sensitive to the criterion chosen. Consequently, in pursuing the goal of 
developing a computationally efficient alternative to solution of the Navier-Stokes 
equations, the potential flow for the geometry in question is presently being used 
as the basic flow to be corrected. Methods for improving or correcting the pressure 
gradients are deferred for a subsequent investigation. 

The implications of this type of pressure approximation may be clarified by 
examination of the inviscid momentum equation written as 


Vp -r 




where q^ = U • U, and vorticity = V x U. In an orthogonal coordinate system, the 

?,-[vp + pV(qV2)j " (12) 


i^ component of Eq . (11) is 
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for small transverse vorticity ^ 2 ^ ^3 small transverse velocities u^, 

(relative to the coordinate system being used), the right-hand side of Eq. (12) is 
of second order and may be neglected. This "small shear" assumption is familiar as 
part of the "secondary flow approximation" (cf. Hawthorne [3]), although it should 
be emphasized that the present analysis accounts for "distortion of Bernoulli 
surfaces" by means of u . 

If the right side of Eq. (12) is neglected, if pressure gradients are derived 

2 - - . 

from the inviscid velocity by setting q = U^. • U^, and if a mean viscous pres- 
sure correction p^(x^) is introduced, the pressure approximation can be written as 

i, • 7p i v[Vp^(X|) - /jV(Uj-up/2] = i, • V(Pj + p^) 


where p^ is the imposed pressure. Typically, for internal flows, p^ is detetmined to 
ensure that an integral mass flux condition is satisfied, such as 



= constant 


(14) 


For external flows, p^ = 0. Arguments favoring this pressure approximation lose some 
validity in strongly curved corner regions, where strong crossflows in the inviscid 
outer region of shear layers are deflected. Nevertheless, the approximation may not 
be seriously in error when considering its overall and somewhat limited role in the 
present context. Finally, this type of pressure approximation is questionable for 
strongly deflected inviscid flows having large transverse vorticity. 

Combining Eqs. (10) and (13) and setting F = F' provides an equation nominally 
governing u^: 

i,-[(U-V)U +(VPy)//5 -V(Uj-Uj)/2] ^3) 

where U is to be written as 


U = + i2(VjU^ + v,) + ijCWjU^ + Wj) 


(16) 


where u^, v^, w^ are components of U^, and v^ , w^ are components of U^. 
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Secondary Vorticity 


An equation governing in compressible flow can be obtained from an approxi- 
mate application of secondary flow theory. The equation governing the growth of 
vorticity along a streamline for incompressible flow with constant viscosity is 
given by Lakshminarayana and Horlock [6] and may be written as 


a 

Js 


( 



^qR 



s • V X Vp + —2 S 

p 


G 


(17) 


where q is velocity magnitude, s is distance along a streamline, ^ is vorticity in 

n ^ 

the direction of the unit principal normal vector n of the streamline, for which R is 
the principal radius of curvature. These quantities are related by the Frenet for- 

A ^ A- 

mula n/R = 9s/ s • Vs. In Eq. (17), G = V x F, and the term containing p vanishes 
if p is constant. Since an intrinsic coordinate system formulation as in Eq. (17) 
provides an "inverse” coordinate system if used to compute U, and since intrinsic 
coordinates are degenerate on no-slip surfaces and nonorthogonal for general rota- 
tional flows, intrinsic coordinates are not attractive for numerical computation. 
However, if the coordinate system used for computation is approximately aligned with 
the flow direction, then an approximate equation governing fi, can be derived from 
Eq. (17) by replacing s by i, as in the following development: 

Xlg s . i| • XI XT, I 

q s • U i|- U U| 

_ ts.Vs).n (i,-Vi,)»n Xin, 

qR s-U iTq ■ U|R, 

where u^ = u^u ; R, is the principal radius of curvature of the x, coordinate, and 
II v l I 

, is vorticity in the direction of n^ , the principal normal of the x, coordinate 

A J- A. J- A A 

line. The quantities n^ and R^^ are defined by the Frenet formula n^^/R^^ = i^^ • Vi^^. 

To illustrate, in an orthogonal coordinate system, 

= - r -^1 (19 

R| L hjhg < 3 Xg axjj 


(18a) 


(18b) 
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where h. , h , denote metric coefficients. If p varies, p is replaced by the 
-L /, o 2 2 ^^ 

imposed pressure p^ as defined in Eq. (13). Finally, taking q ^2. s = i^ in 

the last two terms in Eq. (17), and neglecting streamwise diffusion, Eq. (17) becomes 


u-v 



2X1 


ni 





+ 


I A 


g' 


( 20 ) 


where G’ does not contain streamwise diffusion. 

The transverse vorticity in Eq. (20) contains components which, in orthogonal 

coordinates and assuming an irrotational U^, are given by 


^2 = • VX U 




h, ax, 





( hjWj ) 


(21a) 


du, 




i^.vxu . 


dx - 
2 2 


ax 






(21b) 


In the applications contemplated here, the first term on the right-hand side of each 
of Eqs. (21a-b) is expected to dominate, and the remaining terms may be neglected as 
a convenience, if desired. 

An alternative equation governing vorticity can be derived from the general 
vector equation for vorticity in compressible flow: 


(U- V)X1 - (Xl-V)u + Xl(V-U) + VX — Vp = G 

p 

Streamwise diffusion is neglected in Eq. (22), and if p varies, p is replaced by 
Utilizing continuity and relationships valid for general orthogonal coordinates, 
component of Eq. (22) can be written as 


(22) 

Pi- 

the 


A 

- - 

u- vn, - n - vu, + (Xi,u - U|Xi) • 


a, _ * I 

U-Vp + i -Vx — Vp_ 

P ^ ' P ^ 


* - / 

I • G 


(23) 


15 



The last two terms on the left side of Eq. (23) vanish if p is constant. Because of 
insight provided by previous work on secondary flow theory, Eq. (20) was used instead 
of Eq. (23) in the present application to curved flow geometries. 

Energy Equation 

If the flow being considered is nonadiabatic, solution of an equation governing 
energy is required. The energy equation can be written in a variety of forms, one of 
which is 

pU-VE = V • kVT + U * F + <I> ( 24 ) 

where E is total enthalpy, T is temperature, k is thermal conductivity, and ^ is the 
dissipation function. Solution of Eq . (24) by forward marching integration requires 
only that terms representing streamwise conduction of heat and also streamwise viscous 
diffusion in F be neglected. 

Compressibility Relations 


The foregoing analysis can be applied to incompressible flows simply by setting 
p = p^. Compressibility effects are represented by introducing the perfect gas 
equation of state p = pRT in the imposition of streamwise pressure gradients. For 
moderate subsonic Mach numbers , inviscid pressure gradients can be obtained either 
from a compressible potential flow calculation or from an incompressible potential 
flow corrected for compressibility using either the Prandtl-Glauert formula, Laitone's 
variant [22], or the recent suggestion of Lieblein and Stockman [23]. Replacing p 
in Eq. (13) by the state equation and eliminating temperature using the temperature- 
enthalpy relation 


E 




u-u 

2 


(25) 


where c^ denotes specific heat, the following auxiliary equation relating the imposed 
pressure gradients with density, velocity, and total enthalpy is obtained: 


+ Pv '"l’] 


A 



/ 



(26) 


where y is specific heat ratio. A slight simplification results if U • U is replaced 
by (i^ • U)^ in Eq. (26). 
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In many problems of interest, it can be assumed that the total enthalpy is a 
constant E^. This assumption is reasonable for inviscld flow regions with no heat 
addition and for boundary layers on adiabatic walls provided the Prandtl number is 
unity. To make this assumption, E is replaced by E^ in Eq. (26) , and it is then 
unnecessary to solve the energy equation, even though the flow is compressible. 

Typically, the treatment of nonconstant total enthalpy is of interest for pre- 
dicting heat transfer levels in wall boundary layers. However, it should be noted 
that the geometrical pattern of streamlines for steady inviscid flow of a perfect 
gas depends only on the total pressure gradient and not on the distribution of total 
enthalpy (cf. Hawthorne [3]). Thus, within an assumption of small transverse total 
pressure gradient in the free stream or core region, it is possible to consider core 
flows with two separate streams having large differences in total enthalpy and veloc- 
ity and hence separated by thermal and viscous shear layers, provided the two streams 
have about the same static and total pressure. This observation was used to advantage 
by Kreskovsky, Briley and McDonald [15] in an application of the present analysis to 
flow in turbofan exhaust lobe mixers. 

Viscous and Heat Conduction Terms 


A variety of techniques for suppressing streamwise diffusion is possible, and 
the particular choice in any given application may depend on the complexity one is 
willing to tolerate in the viscous terms. For example, the complete set of viscous 
terms can be written out and all derivatives in the streamwise direction neglected. 
Alternatively, only second derivatives in the streamwise direction may be neglected. 
Here, attention is restricted to incompressible flow but with variable viscosity. In 
the momentum equation (10), F = -V x vV x U , where v is kinematic viscosity. A par- 
ticularly simple expression which neglects contributions from is given by 


A — 

i,* F * 


- I, • Vxi/Vx i|U jUy 


f' 


(27) 


Alternatively, contributions from may be retained as in the approximation 


f, • F = - i , ■ V X X u = F 


(28) 
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Similarly, in the vorticity equation (17), G is defined by G = -V x V x vfi and ap- 
proximated by 


-i|-VxVxi|r/n|»G 


(29) 


Streamwise heat conduction in Eq. (24) is suppressed by the approximation 


V • kVT ^ V* kVgT 


(30) 


Governing System of Equations 


A complete system of six coupled equations governing u^, (|), i/;, E, and p is 

given by Eqs. (8), (9), (15), (24), (26), and (20) or (23). Ancillary relations are 
given by Eq. (5) for composite velocity, Eq. (14) for mass flux, and Eqs. (21a-b) for 
transverse vorticity, Eq. (26) may be omitted if p is constant, and Eq . (24) may be 
omitted for constant E. In the Appendix, these equations are given in general orthog- 
onal coordinates. 
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Numerical Method 


Since techniques for obtaining the basic potential flow solution are well known 
and numerous, they need not be enumerated or discussed here. Instead, the present 
development concentrates on describing the numerical method used to solve the system 
of correction equations. Streamwise derivative terms in the governing equations have 
a form such as Uj^3( )/3x^, and because the streamwise velocity u^^ is very small in the 
viscous dominated region near no-slip walls, it is essential to use implicit algorithms 
which are not subject to stringent stability restrictions unrelated to accuracy 
requirements. Although it is possible to devise algorithms for solution of the cor- 
rection equations as a fully coupled implicit system, such algorithms would require 
considerable iteration for the system of equations treated here, and this would de- 
tract from the overall efficiency. The present method is semi-implicit and seeks to 
reduce the amount of iteration required and yet avoid the more severe stability re- 
strictions of explicit algorithms. The method partitions the system of correction 
equations into subsystems which govern the primary flow, the secondary flow, and the 
turbulence model. The primary-flow subset of equations contains the streamwise 
momentum and energy equations, the state equation, and the integral mass flux rela- 
tion. The secondary- flow subset of equations contains the secondary vorticity 
equation and the scalar and vector potential equations. The turbulence model in the 
present study includes a turbulence kinetic energy equation and a length scale 
relationship. These subsystems are decoupled using an hoc linearization in which 
secondary velocity components and turbulent viscosity are lagged, and are solved 
sequentially during each axial step. 

Summary of Algorithm 

The correction equations are replaced by finite-difference approximations. 
Three-point central difference formulas are used for all transverse spatial deriva- 
tives. An analytical coordinate transformation devised by Roberts [24] is employed 
as a means of introducing a nonuniform grid in each transverse coordinate direction, 
as appropriate, to concentrate grid points in the wall boundary layer regions. 

Second— order accuracy for the transverse directions is rigorously maintained. Two- 
point backward difference approximations are used for streamwise derivatives, although 
this is not essential. 
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To solve the primary flow subsystem of viscous correction equations, scalar ADI 
schemes are used for the momentum and energy equations, and the standard secant itera- 
tion process [25] is used to find the mean pressure drop consistent with the integral 

mass flux relationship and state equation. It is convenient to choose u^ = u^ u as 

1 I V 

the dependent variable for the primary flow momentum equation. 

Given the solution for the primary flow, the secondary flow subsystem can be 
solved. First the scalar potential equation (continuity) is solved using a scalar 
iterative ADI scheme. Next, the secondary vorticity and vector potential equations 
are written as a fully implicit coupled system and solved using an iterative linearized 
block implicit (LBI) scheme (cf. Briley and McDonald [7]). In selecting boundary 
conditions for the secondary flow subsystem, care must be taken to ensure that the 
final secondary velocity satisfies the no-slip condition accurately. Zero normal 
derivatives of tp are specified in the scalar potential equation, and this boundary 
condition corresponds to zero normal velocity. It is not possible to simultaneously 
specify the tangential velocity, however, and thus the 4>-contribution to the secondary 
velocity will have a nonzero tangential (slip) component, denoted v^ , at solid 
boundaries. In the coupled vorticity-vector-potential equations, both normal and 
tangential velocity components can be specified as boundary conditions, since these 
equations are solved as a coupled system. By choosing (a) zero normal velocity and 
(b) -v^ as the ijy-contribution to the tangential velocity, the slip velocity v^ 
arising from the (}) calculation is cancelled, and the composite secondary flow 
velocity including both <f> and t|> contributions will satisfy the no-slip condition 
exactly. 

Finally, the turbulence model requires solution of the turbulence kinetic energy 
equation, and this is done using a scalar ADI scheme. 

A summary of the overall algorithm used to advance the solution a single axial 
step follows. It is assumed that the solution is known at the n level x^ and is 
desired at 

(1) The imposed streamwise pressure gradient distribution is determined 
from an a^ priori inviscid potential flow. 

(2) A value for the mean viscous pressure drop is assumed. Initially, the 

value from the previous step p^ is used. 

n+l 

(3) The momentum equation is solved to determine u^ , and if this is the 
first secant iteration for the mean pressure drop, the energy 
equation is also solved for 
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I 


(4) The density is determined from the equation of state, using the 
imposed pressure with mean correction added, and using and 


(5) Since for internal flows, the integral mass flux relation (14) 

will generally not be satisfied, return to step (2) and repeat 

this process iteratively using the standard secant method [25] 

to find the value of which leads to and satisfying the 

V 1 ° 

integral mass flux relation (14) . The secant method was found in 
practice to converge to five figures on the third iteration. 

ri'4“X 

(6) Using values now available for p and u^^ , the scalar potential 

equation (8) is solved using an iterative scalar ADI scheme, to 
obtain This ensures that the continuity equation is 

satisfied . 


(7) The equations for vorticity (20) and vector potential (9) form 

a coupled system for arid which is solved as a coupled 

system using an iterative LBI scheme. 

(8) Values for the transverse velocities v^ and w^ are computed 
from Eq . (3) . 

(9) The turbulence kinetic energy equation is solved using a scalar 
ADI scheme. Updated values of are then computed from 

and the specified length scale. The turbulence model used here is 
given in the Appendix. 
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axial steps. In Fig. 6, computed results for a typical fully developed radial veloc- 
ity profile are compared with the calculated results of Ghia and Sokhey [12], which 
were made under slightly different similarity assumptions for fully developed flow. 

For fully developed flow, the only difference in the theories are found in the equa- 
tions governing the secondary flows. The present analysis utilizes a streamwise 
vorticity equation derived from viscous secondary flow theory, whereas the secondary 
flow in the results of Ghia and Sokhey is governed by components of the momentum equa- 
tion in the directions of a polar-cylindrical coordinate system. Nevertheless, Fig. 6 
shows that the two predictions are in reasonable agreement. It should be recalled 
that the present analysis permits optional use of the streamwise component of the 
vector vorticity equation in place of the vorticity equation derived from secondary 
flow theory. Use of this equation would presumably reproduce the fully-developed 
results of Ghia and Sokhey. However, both approaches involve coordinate-related 
approximations, and it is not clear at this point which option would provide better 
flow predictions for developing flows. The differences observed in Fig, 6 are small 
and may be attributable in large part to differences in numerical truncation error, 
since different mesh distributions and grid points were used. In Fig. 7, an addi- 
tional comparison is shown for the fully-developed primary flow velocity profile at 
the symmetry line midway between the two endwalls . Again, only minor differences 
between the present solution and that of Ghia and Sokhey are present. These dif- 
ferences presumably reflect the aforementioned differences in secondary flow velocity. 
Also shown in Fig. 7 are the experimental measurements of Mori, Uchida and Ukon [29], 
which do not agree well with either of the computed predictions. The disagreement 
between these measurements and various calculations for fully developed flow has not 
as yet been explained. Figures 8 and 9 contain the present predictions for primary 
and secondary velocity in the fully developed flow region. 

Turbulent Duct Flow Comparison 


The final test case consists of developing turbulent flow in the curved section 
of a 90° bend with straight extensions at the inlet and exit. This case was computed 
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for comparison with the recent experimental measurements of Taylor, \^itelaw, and 

Yianneskis [30] and has a curvature ratio R/H = 2.3, a Reynolds number Re of 40,000 

n 

and a Dean number of 26,370. The potential flow is not a free vortex for this 
geometry, because of the straight sections at the inlet and exit, although a free- 
vortex condition does occur near the center of the bend. For convenience in the 
present calculations, the potential flow pressure gradients were approximated by a 
streamwise distribution of free vortex flows with a variable centerline radius which 
is chosen such that the radial pressure difference between inner and outer radii of 
the bend matches that of an available potential flow calculation for this geometry 
(Humphrey, Taylor and Whitelaw [31]), and it is evident that the present approximation 
closely matches the transverse variation in streamwise pressure gradient present in 
the computed potential flow for this geometry. 

The present calculation was initiated at the start of the bend with boundary 
layer thickness 6 = 0.25H, chosen to match the experimental flow. Because the flow 
begins to turn upstream of the start of the bend, there is secondary flow present at 
the start of the calculation. This secondary flow was estimated from the Squire- 
Winter theory and corrected for no-slip conditions and continuity, as in the previous 
calculation at the Horlock, et al laminar flow. The transverse grid was locally 
refined near walls and consisted of 30 radial points and 20 points for the region 
between the endwall and symmetry plane. To reach the 77.5 degree position at which 
the last measurements in the curved duct were taken, 25 equally spaced axial steps 
were computed, and this required about 4 minutes of CDC 7600 run time. 

Figures 11, 12 and 13 show comparisons at the 30°, 60° and 77.5° positions, 
respectively, of the predicted primary and secondary flow velocities and the measure- 
ments of Taylor, Whitelaw, and Yianneskis [30]. In general, the computed and measured 
velocity profiles are in very good agreement, except as noted below. The radial 
velocity profiles seem to be predicted extremely well, considering the difficulty 
of predicting strong turbulent secondary flows. It should be noted that the velocity 
profiles nearest both the pressure and suction surfaces (r = 0.1 and 0.9, respectively) 
are located well within the boundary layers on these surfaces, and are thus very 
sensitive to radial position. The only significant disagreement between the com- 
puted and measured velocity profiles occurs in the profiles nearest the suction side 
of the duct at the 60° and particularly the 77.5° locations (Figs. 12 and 13). In 
this region, there is considerable distortion of the primary flow as a result of the 
strong vortical secondary flow. The computed secondary flow structure consists of 
very strong radial flow within the endwall boundary layer and toward the inner wall 
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or suction surface of the duct, with radial velocity about 43 percent of the mean 
primary flow velocity. Upon reaching the suction surface, the secondary flow pro- 
ceeds up the suction surface with peak velocity about 30 percent of the primary flow, 
and then near the symmetry plane is deflected away from the suction surface and back 
into the primary flow, forming a vortex structure. The radial flow near the 
symmetry plane is less concentrated than that near the walls and has a peak velocity 
of about 18 percent of the primary flow. 

Since the primary flow is convected by the secondary flow, one effect of the 
strong secondary flow is to distort the primary flow, particularly near the suction 
surface. In Figs. 12 and 13, evidence of this primary flow distortion near the 
suction surface can be seen in both the present solution and in the experimental 
measurements, although agreement between the two is at best qualitative. The present 
computation predicts less distortion of the primary velocity near the suction surface 
than is present in the experimental measurements, particularly in Fig. 13a after 
77.5 degrees of turning. This discrepancy may be due to approximations made in the 
analysis such as the treatment of streamwise pressure gradients or may represent a 
limitation of the turbulence model. Another possible explanation worthy of mention 
is numerical truncation error due to a large axial step size. The step size used in 
this calculation corresponds to a distance of 0.125H along the duct centerline and 
covers about 3 degrees of turning. Since the numerical algorithm employs a first- 
order backward difference formulation with "lagged” secondary velocities for the 
primary momentum and secondary vorticity equations, it may be expected that the 
computed solution would develop more slowly than the actual solution due to the 
first-order "lagging" truncation error. The computed results in Fig. 13a are 
consistent with this behavior, as was verified by comparing the computed primary 
flow velocity profile at r = 0.9 and 0 = 77.5 degrees (Fig. 13a) with the correspond- 
ing measurements at 9 = 60 degrees (Fig. 12a). If the average free stream velocity 
is adjusted to reflect the difference in local static pressure, then the computed 
velocity profile for 77.5 degrees is in excellent agreement with the measured profile 
for 60 degrees. Further study including mesh refinement tests would be required to 
determine whether truncation error has significantly influenced the present compari- 
sons - it is hoped that this mesh refinement study will be performed in the near 
future. Finally, in Fig. 14, the computed primary velocity profiles in the sjnnmetry 
plane midway between the endwall surfaces are compared with the corresponding measure- 
ments. The level of agreement is the same as that present in Figs. 11-13. 
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SUM>IARY A2JD CONCLUSIONS 


(1) The present study includes further development and evaluation of the primary- 
secondary velocity decomposition approach suggested by Briley and McDonald [1] 
for application to viscous subsonic flow in smoothly curved geometries. 

(2) An improved solution algorithm was developed here which provides for solution 
of the equations governing secondary vorticity and vector potential as a 
coupled system, using an iterative linearized block implicit (LBI) scheme. 

This permits implicit specification of no-slip boundary conditions for the 
secondary flow velocities and removes a previous assumption that flow in 

the near wall region is collateral. The present solution procedure is 
effective even with the extreme local mesh resolution which is necessary 
to resolve near-wall sublayer regions in turbulent flow calculations. 

(3) Laminar flow calculations were made for specialized flow cases which would 
permit comparison with other experimental and analytical results. Secondary 
flows predicted by the present method were in good agreement with experimental 
measurements and a boundary layer similarity solution for a laminar flow in a 
strongly curved duct with thin boundary layers. This is a significant com- 
parison and provides an indication that the present analysis adequately 
represents strong cross flows in thin boundary layers when they occur as a 
separate region of the flow distinct from corner flow regions. 

(A) Although fully developed flows have little relevance in applications of 

interest here, the present prediction of a laminar fully developed flow agreed 
well with other fully developed solutions. This provides some validation of 
the approximations used in deriving the equation for secondary vorticity as an 
approximation of secondary flow theory. 

(5) Predictions for turbulent flow in a strongly curved duct with moderate boundary 
layer thickness were compared with measured values of primary and secondary 
velocity. These flow conditions are typical of those occurring in turbo- 
machinery applications of interest here. Considering the complexity of 
this flow, the agreement between predicted and measured velocity is very good 
except near the suction surface after considerable turning. In this region 
there is a strong corner vortex which results in considerable distortion of the 
primary flow. Although the corner vortex and primary flow distortion were 
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predicted in qualitative terms, further evaluation is necessary to 
determine the reason or reasons for the quantitative disagreement in 
this region. 

(6) Overall, the present evaluation suggests that the initial value correction 
method used here does provide a reasonable description of the development 
of thin shear layers and secondary flows in strongly curved ducts. The 
method appears very promising as an economical procedure for making 
detailed predictions for this category of flow problems. 
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APPENDIX A: 


governing equations in orthogonal coordinates 

Velocity Decomposition 

U = i,UjUy + igtVjUy + Vg) + i3(w^u^ + w^) 
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Momentum 
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The transverse vorticity components are given by 
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Pressure Gradient Relation 


Eq. (26) is given by 


where 
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Energy Equation 

To simplify the energy equation, assumptions are made which are appropriate for 
boundary layer flow on walls aligned with either an surface or an - x^ 

surface. Eq . (24) can then be approximated by 
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where the Prandtl number Pr is defined by Pr 
Prandtl number. 


c y/k, and Pr is the turbulent 
P t 
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APPENDIX B: 


TURBUT.ENCE MODEL 


In turbulent flow applications, the stress terms F' and G' in Eqs. (15), (20), 

and (22) contain turbulent shear stresses in the form of fluctuating velocity and 

velocity-temperature correlations. In the present application to flow in curved 

ducts, these turbulent stresses are modeled by introduction of an effective turbulent 

eddy viscosity whose distribution is determined by solution of turbulence-model 

equations. The present treatment employs the one-equation turbulence model developed 

by Shararoth and Gibeling [32]. This model is based upon the turbulence energy partial 

differential equation, an algebraic length scale equation and a relation between dis- 

3 / A 3/2 

sipation, turbulence energy, and the length scale of the form e = k / Jl 

where is a turbulence structural coefficient. To calculate low Reynolds number or 
transitional flows, C is made a function of the turbulent Reynolds number. The 

y 

length scale is determined from a mixing length distribution appropriate for the 
flow under consideration. This turbulence model is related to two-equation models as 
discussed by Launder and Spalding [33], since both models solve a form of the tur- 
bulence kinetic energy equation. The major difference is that the one-equation model 
replaces assumptions necessary to solve the turbulence dissipation equation with 
assumptions regarding the length scale formulation. In the present application, the 
length scale is determined from solution of a momentum integral equation which serves 
to model the strearawise development of a freestream length scale, and from a prescribed 
length scale distribution in cross-sectional planes. 

Using Cartesian tensor notation and overbars to indicate time averaging the 
turbulence kinetic energy is defined as 


k- 


_\_ 

2 ^ 


(B-1) 


and the turbulence dissipation is 


<3 u i (3 u'j 
(3x j 
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The low Reynolds number form of the transport equation for the turbulence kinetic 
energy is given [33] as 


/oU 'Vk= V (^+/z^/cT(^) Vk + p - 2ju. Vk‘/^ • 


(B-3) 


where P is the turbulent production. The constant, is taken as 1.0 as recommended 
by Launder and Spalding [33]. By hypotheses, Prandtl and Kolmogorov has suggested 
that the effective viscosity is proportional to the local density and the product of 
a characteristic turbulent velocity and length scale. The characteristic velocity 
is assumed equal to the square root of the local value of turbulence kinetic energy, 
thus 


= C/oi 


(B-4) 


where c is a constant of proportionality. Through dimensional arguments, the length 
scale may be related to the turbulent kinetic energy and dissipation as 


X- 
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Based on an examination of a large amount of experimental data Launder [33] has 
recommended that c*^ = 0.09 = c^, and under these conditions the length scale may be 
thought of as the conventional mixing length. Using Eq. (5) the turbulence kinetic 
energy equation may be expressed as 

/>U-Vk = V (;i + ^^/o-^) Vk + P - 


For near wall and low Reynolds number flows, c^ is modified by a damping factor. 
Following Shamroth and Gibeling [32] , a turbulence function is defined such that 


4a 
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where a^ is a function of a turbulent Reynolds number originally derived by McDonald 
and Fish [34] for transitional boundary layer flows. 
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In Eq. (B-8), a^ is taken as 0.0115 R^ is the turbulent Reynolds number, and the 
function f(R^ is given as [34] 
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f(R^) = 100. R^ < I 

f(R.r^ " R^^40 


(B-9) 


Between the limits of 1<R^<40 a cubic is used to join the two functional forms given 
by Eq. (B-9). In [34], the turbulent Reynolds number, R^, was defined as an integral 
average across the boundary layer. Here, the low Reynolds number variation of a^^ 
relates to the near wall region rather than to transition from laminar to turbulent 
flow, and R is thus taken as the local ratio of turbulent to laminar viscosity. 
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Use of the one-equation turbulence model requires specification of a length 
scale distribution appropriate for the problem under consideration. The duct flows 
of interest here are presumed to have moderately thin shear layers on the boundary 
walls, and the length scale distribution is thus adapted from previous turbulence 
models for turbulent boundary layers. The mixing length distribution of McDonald and 
Fish [34] is one which has proven effective for a wide range of two-dimensional 
turbulent boundary layers and is easily adapted for present use. This distribution is 
given by 


J. = ( K d/^^) 


(B-11) 


where Z is mixing length, is a free stream mixing length, d is distance from the 

wall, k is the von Karman constant (taken as 0.43), and ^ is a sublayer damping 
function given by 


[ (d"" - 23)/e] 


(B-12) 


+ + 1/2 

Here, P is the normal probability function and d is defined by d = d (x/p) /v, 

where t is shear stress. For equilibrium boundary layers, Jl^ is about 0.09 5, where 
6 is the boundary layer thickness. 

The length distribution of Eq. (B-11) is adapted for present use by taking d as 
distance to the nearest wall and by assigning a representative value at each 
streamwise location in the duct, taken as 0.09 6^ where 6^ is a representative boundary 
layer thickness. The only remaining difficulty is to determine a suitable value of 
Since shear layers occur on all surfaces and in corner regions, and since the 
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free stream primary flow is subject to considerable distortion, it does not appear 
feasible to determine a mean boundary layer thickness from the cross-sectional 
variation of the primary flow velocity. Instead, at the starting location, an 
"effective” boundary layer thickness is defined as 
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li/i 
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where 6^ is the initial boundary layer thickness on "i"th wall, and is the 
length of the "i"th wall of the duct. The boundary layer thickness is then assumed 
to grow in accordance with the momentum integral equation for a two-dimensional 
boundary layer. Assuming a 1/7 power law velocity profile and a pipe flow skin 
friction law consistent with the 1/7 power velocity profile (cf. Schlichting [35]), 
the momentum Integral equation can be written as 
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The free stream velocity gradient du /dx in Eq. (B-14) is determined from the mean 

pressure drop predicted as the forward-marching calculation proceeds. Eq. (B-14) is 

then solved to obtain 6^, and the free stream mixing length is taken as = 0.09 

5 . With the length distribution then available from Eq . (B-11) , the turbulence 

e 

kinetic energy equation (B-6) is solved, and the turbulent viscosity is determined 
from Eq . (B-4) . 
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Blasius Transverse Coordinate, n = y (U/vx) 

Fig. 2 - Yaw Angle Profile for a Circular Arc Duct of 
Rectangular Cross Section, R/H=3.5, X=0.75 H. 
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a. Primary Flow Velocity Contours 



b. Secondary Flow Velocity Vectors 


Fig. 3 - Computed Results for Circular Arc Duct, Re=10^, 
R/H=3.5 at Measuring Station of Horlock et al. 
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from Endwall 


42 



r 

i 


1 


1 

! 

( 


I 

r 
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Fig. 5 - Flow in Surfaces One Grid Point Away from Sidewalls 
("unwrapped*’ to lie in a plane) . 
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Primary Flow Velocity, u/U 



Radial Distance at S3nranetry Line, r 

- Fully Developed Primary Flow Velocity Profile for a 
Circular Arc Duct of Square Cross Section, R/H=14, 
Re=205, K=55. 
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.9 Primary Flow Velocity, u/U 

Fig. 12a - Primary Flow Velocity Profiles After 60 Degrees of Turning in 
a Circular Arc Duct of Square Crossection, R/H=2.3, Re=40000, 
k=26370. 





Primary Flow Velocity, u/U 

Fig. 13a - Primary Flow Velocity Profiles After 77.5 Degrees of Turning in a Circular 
Arc Duct of Square Crossection, R/H=2.3, Re=40000, K=26370. 
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Prediction 

O Measurements of Taylor, Whitelaw, and Yianneskis 


Pressure Side Suction Side 

r = (r-r^)/(r^-r^) = .1 .3 .5 .7 .9 



Radial Plow Velocity, v/U 

Fig. 13b - Radial Flow Velocity Profiles After 77.5 Degrees of Turning in a Circular Arc 
Duct of Square Crossection, R/H=2.3, Re=40000, K=26370. 
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